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The chief use of this table of the zeros and weight factors of the Hermite polynomials 
is in the calculation of integrals over the interval [-co, oo ] ? when the integrand is either 
the product of e~ x2 and a polynomial, or may be closely approximated by e~ x2 times a poly- 
nomial. The zeros and weight factors, x^ n) and a$ n) , respectively, together with the aux- 
iliary quantities $ n) = ai n) exp [(x\ n) ) 2 ], which are useful in computation, are all tabulated 
here for the first twenty Hermite polynomials. The zeros z; w) and /3i n) are given to 13 or 
more decimals, and the weight factors «i U) to 13 significant figures. Although other shorter 
tables have appeared, this present table will enable one to cope with problems requiring 
much higher degree and accuracy, both in problems involving direct quadratures and those 
arising in the numerical solution of linear integral equations for the range [— » , oo ]. Thus 
the use of this table in any direct quadrature can give exact accuracy as far as the 39th 
degree (the only inexactitude arising from the use of rounded values of x^ and aj n) ). 



1. Introduction 

The main purpose of the table of the zeros x^ 
and weight factors (" Christoff el numbers") a$ n) of 
the Hermite polynomials is to provide the basic 
quantities that occur in the quadrature formula 



where 



/: 



Rn- 



f(x)dx = j:a^f(x? ) ) + R n , 



t»/ (2w) (8 

~2 n (2n){2n-l) . . . (n+2)(n+l) 



(1) 



for some £, - «><£< oo [2, p. 101-102, 369]. 2 

When the integrand does not contain e~ x explicitly, 
but behaves like e~* 2 times a polynomial, instead of 
(1), there is 



/»oo /» o 

F (x) dx= 



i=l 

(2) 



For dealing with quadratures of the type in (2), it 
is convenient to follow the example of A. Reiz [1] 
and to tabulate also the quantities p^^af* exp 
[(a^ n) ) 2 L From the form of the remainder term R n , 
(1) is seen to be exact (save for the use of rounded 
values of x^ and c*i n) ) when/(#) is a polynomial 
of at most the (2n— l)th degree. 

2. Numerical Solution of Integral Equations 

Besides problems involving direct quadratures, 
there are those arising in the numerical solution of 
linear integral equations of the second kind for 
integrals over the range [— °°, °°], namely, 



4>(y)=f(y) + *£j(x)K(y,x)dz. 



(3) 



i This project was sponsored in part by the Office of Air Research, AMC, 
USAF. 
2 Figures in brackets indicate the literature references at the end of this paper. 



By considering e x f(x)R(y,x) as a polynomial in x, 
and making use of (2) for y=xf\ i=l, 2, . . . , n, 
the approximate solution of (3) is reduced to the 
solution of a set of only n linear equations for the 
(approximate values of) /(x^). These values of 
/(£- n) ) can then be used in (2) to find/(?/) for any 
value of y. The advantage in the use of (2) instead 
of an equally spaced quadrature formula for solving 
(3) is that here only n points are needed to give 
accuracy obtainable by approximating e x j(x)R(y,x) 
as a polynomial of the (2n— l)th degree in x. For a 
full description of the method, including examples^ 
see A. Reiz [1], especially p. 4 to 10, 16 to 21. 



3. Important Properties of Hermite 
Polynomials 

The Hermite polynomials may be defined by 
the formula 



or by the generating function 






They are given explicitly by the formula 



[n/2] (—lY n \ 



(4) 



(5) 



(6) 



The numerical values of the coefficients of H n (x) for 
the first 30 values of n are contained in H. E. Saizer 
[3]. The polynomials H n (x) satisfy the differential 
equation 



H n " (x) -2xH n '(x) +2nH n (x) =0, 



(7) 
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the recurrence formula 

H n+1 (x)-2xH n (x)+2nH n _ l (x)=0, n^2, (8) 

and the differential-recurrence formula 

H n '(x)=2nH n _ l {x). (9) 

The Hermite polynomials have the orthogonality 
property 



/: 



e- xi H n {x)H m (x)dx=^2 n n\8^. 



(10) 



A most important consequence of (10) is (1), where 
the af are defined by 



a?': 



1 [' e-*H n (x) 

~H n \x?))-* x-xr 



The af are also expressible in the following form 
[2, p. 344], [4, p. 311]: 



a$>> 



wi2 n+1 n\ 



\H n '(xr)) 2 ' 
which, from (9) , can be written as 
t*2— '(n— 1)1 



a? 



(12) 



(13) 



the latter form being more convenient for compu- 
tation. 

4. Description of Tables 

This table gives the zeros x} n) , the weight factors 
a/ n) and ft (n) , for the first 20 Hermite polynomials. 
The absolute values of the zeros, |#/ n) |, are tabulated 
to 15 decimals for n= 1(1)13, to 14 decimals for 
7^=14(1)16, and to 13 decimals for n=17(l)20. To 
each nonvanishing \x£ n) -\ there corresponds the two 
zeros xt ) = + \xi n) \ and x} n) = — \xf n) \. The weight 
factors <4 n) = <x}?li+i are tabulated to 13 significant 
figures. For a given n, zeros having the same ab- 
solute value are associated with the same weight 
factor. The auxiliary quantities ft (n) are tabulated to 
13 decimals. 

5. Method of Computation and Checking 

Since H n (x) is an odd or even polynomial for n odd 
or even, respectively, it was convenient to change the 
variable to X=4x 2 and to calculate Xj n) =4xt )2 , 
where the nonvanishing X^ n) are the zeros of poly- 
nomials A [n i2 ] {X) ) n odd, and B [n/2] (X), n even, of 
degree [n/2] in X. Thus the degree of the polyno- 
mials A [n/2] (X) and B [n/2] (X) did not exceed 10, and 
the coefficients were considerably smaller than those 
of H n (x). A first approximation to Xl n) was ob- 
tained from E. R. Smith [5], which gives the zeros of 



9 d n (e- x ' I/2 ) 

i(x) = (-l)Y 2 / 2 ^ 

by taking twice the square of the latter. The poly- 
nomials A [n/2] {X) and B [n/2] (X) were tabulated for 
five or six values of X in the neighborhood of each 
Xf n) , at intervals of 0.0001 in X The zeros X} n) 
were calculated by a formula for inverse interpolation 
given in H. E. Salzer [6], and then were checked by 
direct interpolation in the polynomials A [n/2] (X) and 
B [n/2] (X), using the Gregory-Newton formula with 
advancing differences. The accuracy of the zero 
Xf n) was determined from the approximate relation: 



. V(n) 0.0001 jp /V(n) . 
error in Xi n) ^ — - — F(X} n) ), 



where F{X} n) ) denotes either A [n/2] (Xl n) ) or 
B[n/2](Xl n) ), and A is the difference between two 
consecutive entries of F(X) on either side of Xj n) . 
The values of x} n) were found from x} n) = 

±— -yJX} n) , and checked back by squaring to get 

Xf n) =4x} n)2 . Additional checks upon the x^ n) were 



performed by calculating 



2>/* } = 



2_n(n— 1) 



i=l 



=(_!)(»- 



■D/2. 



2 

(n+l)l 
((»+l)/2)! 



(14) 
for n odd, (15a) 



where 71-' denotes the omission of the factor x£ n) =0, 
and 



T X} n 
i=l 



■(-!)' 



n 



(n/2) 



for n even. 



(15b) 



The weight factors a} n) were calculated fronT(13), 
using 



H n ^(xi n) )=B [{n . l)/2] (Xn for n odd, 



and 



ffn-iW B) )^V^ T i [ (».i)/2]W ) ) for n even. 
Formula (12) was used to check the a/ n) , in the form 



Oti 



(n)__ 



r*2 tt - 8 n! 



-{x^A\ n/2] (xn} 



2 ; for 7i odd, i^-(n+l) (12a) 



and 



a; 



(n). 



7r^2 n - 3 n! 



-X^{B\ nl2] (XnY 



j for n even. (12b) 



Additional over-all checks upon a/ w) were performed, 
using (1) for P(x) = l, % 2 , A . . . , x 2tn , 2™<2ft,-l< 
2 m+1 , and in a few cases for P(x)=x 2n + 2P , 2 m +2?< 
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2n — l<2 m+1 , the left member of (1) equal to r(&+-j 

for P(x)=x 2k . The choice of these different powers 
of x for P(x) was to ensure that every significant 
figure of the a/ n) was covered in the checks. 

The auxiliary quanitites ft (n) were computed by 

multiplying the weight factors al n) by e x * n /4 , where 

the e x i n Ji were obtained from the NBS tables of the 
exponential function [7] and were checked by the 
relation log ft (n) =log <4 n) +X} n) /4, where the loga- 
rithms were obtained from the NBS tables [8]. 

The functions xf n) , a/ n) , and ft (n) were calculated 
to about two more places than are tabulated, and 
are available in manuscript form. All entries given 
here, which were obtained by rounding the manu- 
script table, are correct to within about a unit in the 
last place. 

The functions in this present table were compared 
with those in three much shorter similar tables given 
by R. E. Greenwood and J. J. Miller [9], A. Reiz [1], 
and Z. Kopal [10], and the few errors of more than a 
unit in the last place in those three tables were noted 
and reported to the journal "Mathematical Tables 
and Other Aids to Computation.' ' Also, the values 
of X/ n) for n= 16 upon the worksheets were compared 
with those in J. Barkley Rosser [11], with perfect 
agreement. Finally, the zeros in the table of the 
Harvard University Laboratory [28] were checked by 
calculating -yj2\xl n) \ } to give complete agreement. 

6. Illustration 

As an example of the use of (1) for quadratures, 
suppose that one wished to evaluate numerically, say 

to 13 places, the integral e~ t2 J (t)dt, where J (t) 

denotes the Bessel function of order zero. It hap- 



pens that that integral is known to be Y 2 ^t e" 1/8 I ([ f > ». 
where I (t) denotes the modified Bessel function of 
order zero (see [26, p. 394]). As the integrand 
e- t2 J (t) is even, and as <4 n) J (xj n) ) = aZ l l i+1 
Jo(—Xi n) ), from (1), whenever n is even, 



/: 



e- t2 J {t)dt^-La^J {xl n) ), 



(16) 



where the summation covers only the positive values 
of x^ n) . In the present example the error in (16) 
(from the term R n in (1) and the inequality 
\Jo m) 0*0 1 = 1) does not exceed 10~ 13 for n as small as 10. 
For the right member of (16), Jo(xf m ) was computed 
to as many places as were sufficient to obtain 
al n) J (x} n) ) to 13 places. The J fe (10) ) was calcu- 
lated by interpolation in the Bessel function table 
of the Harvard Computation Laboratory [27]. The 
following are the results of these computations (the 
extra fourteenth place in a} 10) J (xf 10) ) was carried to 
avoid the cumulation of rounding errors): 



i 


Jo(x^) 


ai (10) /oW 10) ) 


6 
7 
8 
9 
10 


0. 97081 99870 413 
0. 74887 18398 349 
0. 36515 41522 17 

-0.06452 02154 9 

-0. 37050 3331 


0. 59303 76541 6692 

0. 17983 30434 9660 

0. 01236 93757 8926 

-0. 00008 66923 1312 

-0. 00000 28308 0582 



ZX 10) «W 

i=6 



= 0.78515 05503 338. 



This answer differs by only a unit in the thirteenth 
place from the true value of }i^w e~ 1/s Io(l/S) which 
was found (employing the power series for I (t)) to 
be, to 13 decimals, 0.78515 05503 339. 



7. Table of the Zeros and Weight Factors of the First Twenty Hermite Polynomials 



N B) | 


<*!»> 


# n) 


n=l 




n=l 


71=1 




a ooooo ooooo 


ooooo 


1. 7724 53850 906 


1. 77245 38509 


055 


n = 2 




71 = 2 


71 = 2 




0. 70710 67811 


86548 


0. 88622 69254 528 


1. 46114 11826 


611 


ti = 3 




71 = 3 


71 = 3 




0. ooooo ooooo 

1. 22474 48713 


ooooo 

91589 


1. 1816 35900 604 
0. 29540 89751 509 


1. 18163 59006 
1. 32393 11752 


037 
136 


71 = 4 




71 = 4 


71 = 4 




0. 52464 76232 

1. 65068 01238 


75290 

85785 


0. 80491 40900 055 
0. (1)81312 83544 725 


1. 05996 44828 
1. 24022 58176 


950 
958 


n = 5 




71 = 5 


71 = 5 




o. ooooo ooooo 

0. 95857 24646 
2. 02018 28704 


ooooo 

13819 
56086 


0. 94530 87204 829 

0. 39361 93231 522 

0. (1)19953 24205 905 


0. 94530 87204 

0. 98658 09967 

1. 18148 86255 


829 
514 
360 
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\xl n) \ 


o$° 


$F 




n = 6 


71 = 6 


71 = 6 






0.43607 74119 27617 

1. 33584 90740 13697 

2. 35060 49736 74492 


0. 72462 95952 244 

0. 15706 73203 229 

0. (2)45300 09905 509 


0. 87640 13344 362 

0. 93558 05576 312 

1. 13690 83326 745 






71 = 7 


71 = 7 


71 = 7 






0. 00000 00000 00000 
0. 81628 78828 58965 
1. 67355 16287 67471 
2. 65196 13568 35233 


0. 81026 46175 568 

0. 42560 72526 101 

0. (1)54515 58281 913 

0. (3)97178 12450 995 


0. 81026 46175 568 
0. 82868 73032 836 

0. 89718 46002 252 

1. 10133 07296 103 






7i = 8 


71 = 8 


71 = 8 






0.38118 69902 07322 
1. 15719 37124 46780 

1. 98165 67566 95843 

2. 93063 74202 57244 


0. 66114 70125 582 

0. 20780 23258 149 

0. (1)17077 98300 741 

0. (3)19960 40722 114 


0. 76454 41286 517 
0. 79289 00483 864 

0. 86675 26065 634 

1. 07193 01442 480 






71 = 9 


71 = 9 


71 = 9 


1 




0. 00000 00000 00000 

0. 72355 10187 52838 

1. 46855 32892 16668 

2. 26658 05845 31843 

3. 19099 32017 81528 


0. 72023 52156 061 

0. 43265 15590 026 

0. (1)88474 52739 438 

0. (2)49436 24275 537 

0. (4)39606 97726 326 


0. 72023 52156 061 
0. 73030 24527 451 
0. 76460 81250 946 

0. 84175 27014 787 

1. 04700 35809 767 


; 




7i=10 


71=10 


71=10 






0. 34290 13272 23705 
1. 03661 08297 89514 

1. 75668 36492 99882 

2. 53273 16742 32790 
3.43615 91188 37738 


0. 61086 26337 353 

0. 24013 86110 823 

0. (1)33874 39445 548 

0. (2)13436 45746 781 

0. (5)76404 32855 233 


0. 68708 18539 513 
0. 70329 63231 049 
0. 74144 19319 436 

0. 82066 61264 048 

1. 02545 16913 657 






71 = 11 


71=11 


71=11 


^ J 




0. 00000 00000 00000 
0. 65680 95668 82100 
1. 32655 70844 94933 
2. 02594 80158 25755 

2. 78329 00997 81652 

3. 66847 08465 59583 


0. 65475 92869 146 

0. 42935 97523 561 

0. 11722 78751 677 

0. (1)11911 39544 491 

0. (3)34681 94663 233 

0. (5)14395 60393 714 


0. 65475 92869 146 
0. 66096 04194 410 
0. 68121 18810 667 
0. 72195 36247 284 

0. 80251 68688 510 

1. 00652 67861 724 






n=12 


71=12 


71=12 






0. 31424 03762 54359 

0. 94778 83912 40164 

1. 59768 26351 52605 

2. 27950 70805 01060 

3. 02063 70251 20890 
3. 88972 48978 69782 


0. 57013 52362 625 
0. 26049 23102 642 
0. (1)51607 98561 588 
0. (2)39053 90584 629 
0. (4)85736 87043 588 
0. (6)26585 51684 356 


0. 62930 78743 695 
0. 63962 12320 203 
0. 66266 27732 669 
0. 70522 03661 122 
0. 78664 39394 633 
0. 98969 90470 923 


\ 




71=13 


71 = 13 


71 = 13 






0. 00000 00000 00000 

0. 60576 38791 71060 

1. 22005 50365 90748 

1. 85310 76516 01512 

2. 51973 56856 78238 

3. 24660 89783 72410 
4. 10133 75961 78640 


0.60439 31879 211 
0. 42161 62968 985 
0. 14032 33206 870 
0. (1)20862 77529 617 
0. (2)12074 59992 719 
0. (4)20430 36040 271 
0. (7)48257 31850 073 


0. 60439 31879 211 
0. 60852 95837 033 
0. 62171 60552 868 
0. 64675 94633 158 
0. 69061 80348 378 
0. 77258 08233 517 
0. 97458 03956 399 






71=14 


71=14 


71=14 






0. 29174 55106 7256 
0. 87871 37873 2940 
1.47668 27311 4114 
2. 09518 32585 0772 

2. 74847 07249 8540 

3. 46265 69336 0227 

4. 30444 85704 7363 


0. 53640 59097 121 
0. 27310 56090 642 
0. (1)68505 53422 347 
0. (2)78500 54726 458 
0. (3)35509 26135 519 
0. (5)47164 84355 019 
0. (8)86285 91168 125 


0. 58406 16905 220 
0. 59110 66670 432 
0. 60637 97391 261 
0. 63290 06064 723 
0. 67770 67591 924 
0. 75998 70873 976 
0. 96087 87030 257 


- 


K 
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« 



Wl 


d^ 


#° 


71 = 15 




71 = 15 


71=15 


0. 00000 00000 


0000 


0. 56410 03087 264 


0. 56410 03087 264 


0. 5650G 95832 


5558 


0. 41202 86874 989 


0.56702 11534 466 


1. 13611 55852 


1092 


0. 15848 89157 959 


0. 57619 33502 835 


1. 71999 25751 


8649 


0. (1)30780 03387 255 


0. 59302 74497 642 


2. 32573 24861 


7386 


0. (2)27780 68842 913 


0. 62066 26035 270 


2. 96716 69279 


0560 


0. (3)10000 44412 325 


0. 66616 60051 091 


3. 66995 03734 


0445 


0. (5)10591 15547 711 


0. 74860 73660 169 


4. 49999 07073 


0939 


0. (8)15224 75804 254 


0. 94836 89708 276 


71=16 




71=16 


71=16 


0.27348 10461 


3815 


0. 50792 94790 166 


0. 54737 52050 378 


0. 82295 14491 


4466 


0. 28064 74585 285 


0. 55244 19573 675 


1. 38025 85391 


9888 


0. (1)83810 04139 899 


0. 56321 78290 882 


1. 95178 79909 


1625 


0. (1)12880 31153 551 


0. 58124 72754 009 


2. 54620 21578 


4748 


0. (3)93228 40086 242 


0. 60973 69582 560 


3. 17699 91619 


7996 


0. (4)27118 60092 538 


0. 65575 56728 761 


3. 86944 79048 


6012 


0. (6)23209 80844 865 


0. 73824 562^2 777 


4. 68873 89393 


0582 


0. (9)26548 07474 Oil 


0. 93687 44928 841 


n = 17 




71=17 


71=17 


0. 00000 00000 


000 


0. 53091 79376 249 


0. 53091 79376 249 


0. 53163 30013 


427 


0. 40182 64694 704 


0. 53307 06545 736 


1. 06764 87257 


435 


0. 17264 82976 701 


0. 53976 31139 085 


1. 61292 43142 


212 


0. (1)40920 03414 976 


0. 55177 73530 782 


2. 17350 28266 


666 


0. (2)50673 49957 628 


0. 57073 92941 245 


2. 75776 29157 


039 


0. (3)29864 32866 978 


0. 59989 27326 678 


3.37893 20911 


415 


0. (5)71122 89140 021 


0. 64629 17002 129 


4. 06194 66758 


755 


0. (7)49770 78981 631 


0. 72874 83705 871 


4. 87134 51936 


744 


0. (10)45805 78930 799 


0. 92625 41399 895 


n=18 




71=18 


71=18 


0. 25826 77505 


191 


0. 48349 56947 255 


0. 51684 58364 816 


0. 77668 29192 


674 


0. 28480 72856 700 


0. 52063 49466 761 


1. 30092 08583 


896 


0. (1)97301 74764 132 


0. 52858 94429 188 


1. 83553 16042 


616 


0. (1)18640 04238 754 


0. 54157 86786 621 


2. 38629 90891 


667 


0. (2)18885 22630 268 


0. 56127 90455 498 


2. 96137 75055 


316 


0. (4)91811 26867 929 


0. 59095 30034 631 


3.57376 90684 


863 


0. (5)18106 54481 093 


0. 63763 01720 062 


4.24811 78735 


681 


0. (7)10467 20579 579 


0. 71999 33831 053 


5.04836 40088 


745 


0. (11)78281 99772 116 


0. 91639 35375 519 


71=19 




71=19 


71=19 


0. 00000 00000 


000 


0. 50297 48882 762 


0. 50297 48882 762 


0. 50352 01634 


239 


0. 39160 89886 130 


0. 50461 53313 522 


1. 01036 83871 


343 


0. 18363 27013 070 


0. 50967 93750 983 


1. 52417 06193 


935 


0. (1)50810 38690 905 


0. 51863 31937 003 


2.04923 17098 


506 


0. (2)79888 66777 723 


0. 53240 23605 516 


2.59113 37897 


945 


0. (3)67087 75214 072 


0. 55269 46209 671 


3. 15784 88183 


476 


0. (4)27209 19776 316 


0. 58277 95297 628 


3. 76218 73519 


640 


0. (6)44882 43147 223 


0. 62965 66326 368 


4. 42853 28066 


038 


0. (8)21630 51009 864 


0. 71188 18743 371 


5.22027 16905 


375 


0. (11)13262 97094 499 


0. 90719 87960 928 


71 = 20 




n = 20 


71 = 20 


0. 24534 07083 


009 


0. 46224 36696 006 


0. 49092 15006 667 


0. 73747 37285 


454 


0. 28667 55053 628 


0. 49384 33852 721 


1. 23407 62153 


953 


0. 10901 72060 200 


0. 49992 08713 363 


1. 73853 77121 


166 


0. (1)24810 52088 746 


0. 50967 90271 175 


2. 25497 40020 


893 


0. (2)32437 73342 238 


0. 52408 03509 486 


2. 78880 60584 


281 


0. (3)22833 86360 163 


0. 54485 17423 644 


3. 34785 45673 


832 


0. (5)78025 56478 532 


0. 57526 24428 525 


3. 94476 40401 


156 


0. (6)10860 69370 769 


0. 62227 86961 914 


4. 60368 24495 


507 


0. (9)43993 40992 273 


0. 70433 29611 769 


5. 38748 08900 


112 


0. (12)22293 93645 534 


0. 89859 19614 532 
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